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Abstract 

We study the amplitude of deeply virtual Compton scattering in next-to-leading order of 
perturbation theory including the two-loop evolution effects for different sets of skewed parton 
distributions (SPDs). It turns out that in the minimal subtraction scheme the relative radiative 
corrections are of order 20-50%. We analyze the dependence of our predictions on the choice 
of SPD, that will allow to discriminate between possible models of SPDs from future high pre- 
cision experimental data, and discuss shortly theoretical uncertainties induced by the radiative 
corrections. 
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Virtual Compton scattering (VCS) in a light cone dominated region, — usually referred to as 
the deeply virtual Compton scattering (DVCS) [0-0, — is a favourable process to get access to 
the so-called skewed parton distributions (SPDs) |l|-[f§. The latter can be viewed as a general- 
ization of the conventional parton densities as measured in deep inelastic scattering (DIS). They 
contain complementary information about the internal structure of hadrons, e.g. the total angular 
momentum carried by quarks and gluons ||. 

The study of SPDs promises to become one of the main issues of hadronic physics for the 
next decade. This marks general trends in the development of experimental techniques towards 
exclusive reactions at the facilities like CERN and DESY. In general, it will be a delicate task to 
measure DVCS, since it is contaminated by the Bethe-Heitler (BH) process. However, it is rather 
encouraging that a first DVCS signal has already been detected by the ZEUS collaboration J/]] . In 
addition, azimuthal, spin and charge asymmetries || allow to get access to the interference term 
between DVCS and BH processes and this should allow to separate real and imaginary parts of 
the DVCS amplitude. Such measurements would pin down plausible shapes of SPDs. 

The theoretical goal is, therefore, to push our understanding of SPDs to the level reached 
nowadays for the usual parton distribution functions. One of the most important questions is the 
Q 2 -evolution of SPDs. In earlier papers we gave a complete solution at the NLO level ||. In 
the present paper we study the role of the NLO radiative corrections for the DVCS amplitude 
quantitatively in the region xb^ > 10~ 2 including the effects of scaling violations. Our aim is to 
decide whether an extension to NNLO is necessary to obtain a stable interpretation of a given set 
of DVCS data. For a first rough comparison with experimental measurements NLO is probably 
fine but to judge the potential of SPDs measurement the precise size of the corrections as well as 
the uncertainties involved have to be known. 

The VCS amplitude is defined in terms of the time ordered product of two electromagnetic 
currents 



Here P\ and P 2 are the momenta of the initial and the final hadrons, respectively. The incoming 
photon with momentum q\ has a large virtuality. The scaling variables, £ and rj, which allow to 
describe different "two-photon" processes in the light-cone dominated region are introduced as 
follows 



In DVCS kinematics — q\ ^> m1 admn , q\ = and we have in fact only one scaling variable £ since 



T flv {P,q 1 ,q 2 )=i J dxe ix «{P 2 ,S 2 \T{j fl {x/2)j v {-x/2)}\P 1 ,S 1 ). 
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From the experimental point of view it is more appropriate to work with 
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the variables — q\ and xbj = — g 2 /2(Pigi), which are related to the variables (0) by 

In the kinematical region we are interested in there are two leading twist contributions 

T^P, qi, q 2 ) = -9lfi& V = £, Q\ A 2 ) + >W^ifc r/ = Q 2 , A 2 ) + ■ • • , (4) 

where the ellipsis stand for a leading twist-two contribution coming from longitudinally polarized 
photons which, however, does not appear in the DVCS cross section, as well as higher twist func- 
tions. The transverse part of the metric tensor, denoted by and the e-tensor are contracted 
with the projection operators V a p = g a p — qiaQip I '(Q1Q2) ■ Therefore, current conservation is man- 
ifest. Our definitions are chosen in such a way that in the forward case, i.e. A = 0, the usual 
structure functions measured in DIS are 



Fi(x Bi ,Q 2 ) = ^-lm^ = x Bi ,0,Q 2 ,0), gi{x Bh Q 2 ) = -Um&tf = x Bj , 0, Q 2 , 0). (5) 
Ztt Ztt 



Since the virtuality of the incoming photon is deep in the Euclidean domain the hadron is 
probed almost on the light cone x 2 ~ 0. Hence, the amplitude can be straightforwardly treated 



by a non-local version of the light-cone operator product expansion p0| . More recently, it has 
been proven that assuming a smooth SPD the collinear singularities are indeed factorizable at 
leading twist to all orders of perturbation theory || [II], |I^] . Therefore, the amplitudes T\ and Q\ 
factorize in a perturbative hard scattering amplitude and a SPD q(t, r], A 2 , /j?) as 



(£, v , t, Q 2 , /i 2 ) V(*, V, A 2 , fi 2 ) + ^r- G T a (t 1> *> ^> ^ 1> A '> A* 2 ) 

N fV 



(6) 

where T v = T\ y T A = Q\ and the sum runs over quark species Q = u,d, s with electrical charge 
eg- 

The non-perturbative input is concentrated in the SPDs which read in the Leipzig conventions 

E 

(QqV \(t, V ) = [pe MP +(P 2 S 2 \ij(-Kn){ 7+ WMIPA), (7) 



7+75 



{Xa}^) = ^ / ^e^+<J- 2 ^ 2 |^(-^) |^^ + |G'- + (^)|i- 1 S' 1 >. (8) 

For the quark distributions we have omitted for brevity the flavour indices. A form factor de- 
composition would give us the functions H, E and H, E for the parity even and odd sectors, 
respectively, as introduced in Ref. B. 
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Using the spatial parity and scaling properties evident from the convolution-type formulae of 
Ref. |13| we write the perturbative expansion of T as 

\1] 1] J Z7T \T] J] H Z J V 7 

where the "—"("-)-") sign corresponds to the V(A) channel. The tree level coefficient functions 
are 

Q T (o) (£,1) = ) , G T<°) = 0. (10) 

The hard scattering amplitude can be calculated by making use of standard methods of 
perturbative QCD ||14j| . At the same time there exists an interesting possibility to predict |L3 



these quantities with the help of the conformal operator product expansion (COPE). It has been 
introduced more than two decades ago in Ref. [|T5[] and consequently applied to exclusive processes 
at leading order ||16|| . The main advantage of the COPE is that under the assumption of conformal 
covariance it predicts the Wilson coefficients of local conformal operators with a given conformal 
spin entering the expansion of the product of two currents up to a normalization constant. The 
latter is fixed by the known Wilson coefficients of forward DIS. Since the COPE is valid only for 
|£| > 1 we perform a summation over all conformal partial waves and obtain a representation of 
the hard scattering amplitudes as a convolution of kernels. E.g. for the non-singlet case we have 

T(£,r/,t,QV) = -F ® (%Y ^ ®C (-,-;a s ) , (11) 

where the convolution is defined for any test functions Tx(x,y) and T 2 (x,y) with t% ® r 2 = 
/ rjjrTi(x, y)T 2 (y, z). The evolution kernel V f^, and the coefficient function C (jj, ^'i a sj ar e 
diagonal w.r.t. the conformal waves — the Gegenbauer polynomials: 

and their eigenvalues 7j(a s ) and Cj(a s ) coincide with anomalous dimensions and Wilson coeffi- 
cients, respectively, which appear in DIS. A closed form of the function F can be deduced order 
by order in a s from the sum 

2x 



F(T , A V f 2X Bti + hj + 2) /l+j+7,/2, 2+ J+7j /2 



3=0 



1+X 



Cf(y). (13) 



The expressions (|Tl]-|T^) which are exact in QCD up to NLO0 define the so-called conformal 

subtraction (CS) scheme. To obtain the MS results one has to perform a scheme transformation 

2 Beyond this order the trace anomaly appears and provides conformal noncovariant terms proportional to 
0/g[a s /(2w) + O(a 2 s )}. 
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[O, [L7|, which is governed by a conformal anomaly appearing in special conformal Ward identities 



18| . Finally, the NLO coefficient functions read 
1 
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We have used here the scaling property mentioned above and the fact that the hard scattering 
amplitude for DVCS, i.e. £ = rj, is a simple analytical continuation of the ones for the production 
of a scalar and pseudo-scalar mesons, respectively, in the collision of a real and a highly virtual 
photon. 

To the same accuracy we have to include the two-loop evolution of the SPDs. Based on the 
knowledge of NLO anomalous dimensions constructed from conformal constraints we will use 
the methods of Ref. || where we have already dealt with scaling violations in NLO. 

For numerical estimates of the DVCS amplitude we expand the hard scattering amplitude 
in terms of Legendre polynomials and reexpress the expansion coefficients in terms of conformal 
moments. This allows us to include easily the NLO evolution of the SPDs as described in Ref. 
[|9] . Unfortunately, this method is reliable only for moderate values of XBj since the convergence of 
the series of polynomials is not under control for very low Xb$- Already for xb] ~ 0.05 one needs 
about 180 polynomials. 

To model the SPD one can first explore the simplest possibility of equating it to the usual 
forward parton density. We will designate this as FPD-model. For small rj, where the so-called 
DGLAP region with \t\ > \rj\ dominates, this may be justified. However, it is an open problem 
how the exclusive region \t\ < \rj\ corresponding to the production or absorption of a meson-like 
state looks like. To get a more realistic model one can use the relation of SPDs to the so-called 
double distributions (DDs) f(z-,z + ), introduced in |l| and rediscovered in 



q(t,v,Q 2 



dz 4 



1-1*4-1 



dz^5(z + + r/z^ — t)f(z- 



Q 2 



(18) 



The latter, according to Ref. |19[], is given by the product of a forward distribution f(z) (more 
precisely q(z) for quarks and zg(z) for gluons) with a profile function ir 



(19) 



4 




1.25 



rH|n 

to 



0.75 




0.25 



Figure 1: In (a) and (b) we show the SPD for the sum of the non-polarized u and u quarks and 
the gluon densities, respectively. The FPD-model (dotted line) and the DD-model with MRS pOf 
parametrization (dashed line) are modeled at a scale of 4 GeV 2 , while the DD-Model with GRV 



2l| parametrization (solid line) has been taken at the momentum scale 0.4 GeV and has been 
evolved afterwards with NLO formulae to Qq = 4 GeV 2 . The skewedness parameter is r] = 1/3 
which corresponds to XBj — 1/2. 



where tt for quarks and gluons is given by 

3[1-M] 2 -^ 



7T U-, Z 4 



i5{[i-M 2 -^ 2 }' 



4 [l-\z + \f ' 16 [l-k + |] 5 ' (20) 

respectively. We will refer to this prescription as the DD-model. Note, that parton distributions 
used here have the support — 1 < z < 1 and, therefore, the contributions —q(—z), +Aq(—z) with 
z > have the usual interpretation as antiquarks. 

Now we are in a position to study the NLO amplitudes in the MS scheme for the models 
specified above. In what follows we use the convention Q 2 = —q\ and equate the factorization scale 
with -q\. The value of A^g is A^ G = 246 MeV. Since we will rely on available parametrizations 
of the parton densities which are defined at different scales the DD-models will differ as well. 
Given a forward distribution at an input scale e.g. Qq = 4 GeV 2 , as e.g. for the MRS fit [ 2(| , 
it can be viewed as evolved only according to the DGLAP equation from a lower scale and then 
folded with a profile to form the DD-model. If the former is defined at an input of Qq = 0.4 GeV 2 , 
like e.g. the GRV densities JO], to confront both models we have to evolve GRV-based SPDs with 
non-forward evolution equations. We demonstrate these features of the DD-models in Fig. [I] where 



we have compared the MRSSO |2(J with the GRV parametrization |21J according to the procedure 
sketched above. 

In general we have found that the amplitudes calculated from the models mentioned above 
are similar in shape. However, in the parity even sector the DVCS amplitude is very sensitive to 
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the z — > behaviour of the sea quarks. Thus, different models will produce qualitatively different 
predictions. This is demonstrated in Fig. 0(a) for T\ y where even the sign of the amplitude for 
small XBj differs for the FPD- and DD-model. Note that the sign of the FPD will change for very 
low a;Bj and that a slightly different parametrization gives a prediction similar to the one for a 
DD-model. We see from this figure that the NLO corrections are as large as 50% and even more. 
Evidently, the imaginary part is more sensitive to the shape of the model distribution than the real 
part. Their ratio is shown in Fig. 0(b). In Fig. ||(c) we compare the predictions of the DD-model 
for the MRS and the GRV parametrizations taken at different input scales as explained above. 
This affects the size of both the real and imaginary part, so that the ratio only slightly changes 
with the scale as shown in Fig. 0(d). The ratio is also not sensitive to the radiative corrections 
in the coefficient functions. Note that the evolution will suppress the DGLAP region, so that the 
asymptotic distributions are concentrated in the region \t\ < rj and are given by the terms with 
the lowest conformal spin in the conformal expansion of SPD. Thus for asymptotically large Q 2 
the imaginary part has to vanish and the radiative corrections to the real part are determined by 
the lowest Wilson coefficient in the COPE. 

In Fig. we show the predictions for Q\ with the GSA parametrization |22| and investigate 
the size of radiative corrections in detail. In this case we have found a similar model dependence 
as in the previous case: predictions are sensitive to the sea quark parametrization which turn on 
at XBj ~ 0.1 and also may cause very substantial radiative corrections. Since the polarized sea 
quark distributions are not well known, we simply neglect them. The real and imaginary parts 
for the DD-model are shown in Fig. 0(a,b), respectively, at the input scale 4 GeV 2 and evolved 
upwards to the scale 10 GeV 2 . Not surprisingly the size of the radiative corrections is decreasing 
with increasing Q 2 . Note that for large xbj the ratio of real to imaginary part provides again a 
useful information about the SPD as discussed above. The size of radiative corrections for T\ and 
Qi are similar. In Fig. 0(c,d) we demonstrate the factorization scale dependence, where we took 
in the hard scattering amplitude and SPD the scale to be /i 2 = {<2 2 /2, Q 2 ,2Q 2 }. We have also 
studied the perturbative corrections in the conformal subtraction scheme, however, we have found 
that in the DVCS kinematics this scheme only slightly reduces the NLO corrections. 

To summarize, we have studied the DVCS amplitude in perturbative QCD for the region 
^Bj > 0.05. It turns out that a separate measurement of real and imaginary part in this kinematic 
domain would provide us with information about the SPDs, which is indispensable to discriminate 
between different models. Moreover, the amplitude at xbj ~ 0.1 depends in a crucial way on both 
the model and the used parametrization of the forward parton distributions, especially, of the 
sea distribution. Generally, we found that the radiative correction can be as large as 50% and 
more, depending on XBj- Fortunately, the ratio of real to imaginary part is less sensitive to these 
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Figure 2: The real part (a) and the ratio of the real to imaginary part (b) of J-\, respectively, for 
the DD- and FPD-model with the MRS parametrization taken at the input scale Qq = 4 GeV 2 . 
Here the DD(FPD)-model at LO and NLO is shown as dashed (dotted) and solid (dash-dotted) 
lines, respectively. Plots (c) and (d) show the same as in (a) and (b), respectively but for MRS 
and GRV parametrizations. Here MRS (GRV) amplitude is plotted at LO as dashed (dotted) and 
at NLO as solid (dash-dotted) lines, respectively. 



radiative corrections and is only mildly scheme dependent. Thanks to this property it can give us 
more insights into the structure of SPDs. 

The kinematical situation relevant for HERA experiments requires an extension of our anal- 
ysis downwards to very low xbj- In this case the polynomial method used here will not be the 
most appropriate tool. Fortunately, a direct numerical convolution of the hard scattering part 
and skewed parton distributions can be carried out without major difficulties. The remaining 
problem is to evolve the SPD in this kinematics, which may be solved by numerical integration 
of the evolution equation. Recently, the evolution kernels at two-loop order have been completely 
constructed [2^] from the knowledge of conformal anomalies and splitting functions and extensive 
use of supersymmetric constraints. This provides in future the opportunity to study evolution 
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Figure 3: The real (a) and the imaginary (b) part of Q\ for the DD-model with the GS p2| 
parametrization taken at the input scale Qq = 4 GeV 2 [LO (dotted line) and NLO (solid line)] as 
well as evolved to Q 2 = 10 GeV 2 [LO (dash-dotted line) and NLO (dashed line)]. Plot (c) and (d) 
shows the relative radiative corrections 5 = [(NLO — LO) /LO] again for the real and imaginary 
part, respectively for Q 2 = 4 GeV 2 . The factorization scale jj 2 is set equal to <2 2 /2 (dotted line), 
Q 2 (solid line), and 2Q 2 (dash-dotted line). 



effects also at very low XBj which is a very important task in view of ongoing experiments on 
diffractive meson production at HERA. 

To decrease the theoretical uncertainties due to radiative corrections, it is necessary to go 
beyond NLO. In a first step one can include only the hard-scattering amplitude to two-loop order 
accuracy. Of course, a direct calculation will be very cumbersome. Fortunately, however, in the 
conformal limit of the theory, when the QCD Gell-Mann-Low function formally equals zero, we 
can rely directly on the COPE and it seems feasible to extend this technology to NNLO. Conformal 
symmetry breaking effects can be perturbatively calculated and one piece of information is already 
available ||24[ . From the practical point of view it should be stressed that all high hopes connected 
to SPDs rely on the assumption that the higher order terms are controllable. We want to emphasize 
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that our approach allows to study efficiently the structure of perturbative corrections and might 
be crucial to establish the same rigour of treatment as for the forward parton densities. 
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